Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



New Limits on the Possible Variation of Physical Constants 



M. J. Drinkwater 1 , J. K. Webb 1 , J. D. Barrow 2 , V. V. Flambaum 

1 School of Physics, University of New South Wales, Sydney 2052, Australia 

2 Astronomy Centre, University of Sussex, Brighton, BN1 9QH 



Accepted for publication in MNRAS 



> 

o 
in 



> 
o 

OS 
(N 



Of 

6 

C3 



ABSTRACT 

Recent detections of high-redshift absorption by both atomic hydrogen and molecular 
gas in the radio spectra of quasars have provided a powerful tool to measure possible 
temporal and spatial variations of physical 'constants' in the Universe. 

We have compared the frequency of high-redshift Hydrogen 21 cm absorption with 
that of associated molecular absorption in two quasars to place new (1 sigma) upper 
limits on any variation in y = g p a 2 (where a is the fine structure constant and g p is the 
proton g-factor) of \ Ay/y\ < 5x 10~ 6 at redshifts z = 0.25 and z = 0.68. These quasars 
are separated by a comoving distance of 3000 Mpc (for H = 75kms~ 1 Mpc~ 1 and 
qo = 0). We also derive limits on the time rates of change of \g' p /g p \ < 1 x 10~ 15 y" 1 and 
\a/a\ < 5 x 10 _16 y _1 between the present epoch and z = 0.68. These limits are more 
than an order of magnitude smaller than previous results derived from high-redshift 
measurements . 
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1 INTRODUCTION 

Quasar absorption systems present ideal laboratories in 
which to search for any temporal or spatial variation in 
the assumed fundamental constants of Nature. Such ideas 
date back to the 1930s, with the first constraints from spec- 
troscopy of QSO absorption systems arising in the 1960s. 
An historical summary of the various propositions is given 
in Varshalovich & Potekhin (1995) and further discussion of 
their theoretical consequences is given in Barrow & Tipler 
(1986). 

The subject of varying constants is of particular current 
interest because of the new possibilities opened up by the 
structure of unified theories, like string theory and M-theory, 
which lead us to expect that additional compact dimensions 
of space may exist. Although these theories do not require 
traditional constants to vary, they allow a rigorous descrip- 
tion of any variations to be provided: one which does not 
merely 'write in' the variation of constants into formulae 
derived under the assumption that they do not vary. This 
self-consistency is possible because of the presence of extra 
dimensions of space in these theories. The 'constants' seen 
in a three-dimensional subspace of the theory will vary at 
the same rate as any change occurring in the extra compact 
dimensions. In this way, consistent simultaneous variations 
of different constants can be described and searches for vary- 
ing constants provide a possible observational handle on the 
question of whether extra dimensions exist (Marciano 1984, 
Barrow 1987, Damour & Polyakov 1994). 



Prior to the advent of theories of this sort only the time 
variation of the gravitational constant could be consistently 
described using scalar-tensor gravity theories, of which the 
Brans-Dicke theory is the simplest example. The modeling 
of variations in other 'constants' was invariably carried out 
by assuming that the time variation of a constant quantity, 
like the fine structure constant, could just be written into the 
usual formulae that hold when it is constant. The possibility 
of linked variations in low-energy constants as a result of 
high-energy unification schemes has the added attraction of 
providing a more powerful means of testing those theories 
(Marciano 1984, Kolb, Perry & Walker 1986, Barrow 1987, 
Dixit & Sher 1988, Campbell & Olive 1995). 

Higher-dimensional theories typically give rise to rela- 
tionships of the following sort 

ai(m f ) = AiGn%, = Bi\ n (£ p i / R) m ; n, m constants (1) 

7T _1 C t j[ln(m,/m 3 ) + 6(fi — rrij) ln(raj//«)] + A; 

where cti(..) are the three gauge couplings evaluated at the 
corresponding mass scale; fi is an arbitrary reference mass 
scale, m* is a characteristic mass scale defining the theory 
(for example, the string scale in a heterotic string theory); 
A is some dimensionless string coupling; £ p i = G~ 1//2 is the 
Planck length, and R is a characteristic mean radius of the 
compact extra-dimensional manifold; dj are numbers de- 
fined by the particular theories and the constants Ai and Bi 
depend upon the topology of the additional ( > 3) dimen- 
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sions. The sum is over j — leptons, quarks, gluons, W , Z 
and applies at energies above /j, ~ 1 GeV (Marciano 1984). 
The term Ai corresponds to some collection of string thresh- 
old corrections that arise in particular string theories or an 
over-arching M theory (Antoniadis & Quiros 1997). They 
contain geometrical and topological factors which are spec- 
ified by the choice of theory. By differentiating these two 
expressions with respect to time (or space), it is possible to 
determine the range of self-consistent variations that are al- 
lowed. In general, for a wide range of supersymmetric unified 
theories, the time variation of different low-energy constants 
will be linked by a relationship of the form (where we con- 
sider to denote the time derivative of etc.) 

5(A =<*!—+> S 2i -^+S 3 — + > &y — + <5 6 -r + ...(2) 
3 G ^— ' oq m* ^— ' rnj A 

where = m e /m p is the electron-proton mass ratio. It is 
natural to expect that all the terms involving time deriva- 
tives of 'constants' will appear in this relation unless the 
constant S prefactors vanish because of supersymmetry or 
some other special symmetry of the underlying theory. This 
relation shows that, since we might expect all terms to be 
of similar order (although there may be vanishing constant 
8 prefactors in particular theories), we might expect varia- 
tions in the Newtonian gravitational 'constant', G/G, to be 
of order a /a 2 . 

In this paper, we consider only the bounds that can 
be placed on the variation of the fine structure constant and 
proton g factor from radio observations of atomic and molec- 
ular transitions in high redshift quasars. To do this we ex- 
ploit the recent dramatic increase in quality of spectroscopic 
molecular absorption at radio frequencies, of gas clouds at 
intermediate redshift, seen against background radio- loud 
quasars. Elsewhere, we will consider the implications of si- 
multaneous variations of several 'constants' and show how 
these observational limits can be used to constrain a class 
of inflationary universe theories in which small fluctuations 
in the fine-structure constant are predicted to occur. 

The rotational transition frequencies of diatomic 
molecules such as CO are proportional to h/(Ma 2 ) where M 
is the reduced mass and a — h 2 /(m e e 2 ) is the Bohr radius. 
The 21 cm hyperfine transition in hydrogen has a frequency 
proportional to /x p /is/(?la 3 ), where fx p — g p eh/ '(4m p c), g p is 
the proton g-factor and \ib = eh/(2m e c). Consequently (the 
dependence of the results on the proton-neutron mass differ- 
ence m p —m n is very small so we assume m p /M is constant) 
the ratio of a hyperfine frequency to a molecular rotational 
frequency is proportional to g p a 2 where a = e 2 /(he) is the 
fine structure constant. Any variation in y = g p o? would 
therefore be observed as a difference in the apparent red- 
shifts: Az/(1 + z) ~ Ay/y. Redshifted molecular emission is 
hard to detect but absorption can be detected to quite high 
redshifts (see review by Combes & Wiklind 1996). Recent 
measurements of molecular absorption in some radio sources 
corresponding to known HI 21 cm absorption systems give us 
the necessary combination to measure this ratio at different 
epochs. 

Common molecular and HI 21 cm absorptions in the 
radio source PKS 1413+135 have previously been studied 
by Varshalovich & Potekhin (1996). They reported a dif- 
ference in the redshifts of the CO molecular and HI 21 cm 
atomic absorptions which they interpreted as a mass change 
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021 2+7^ 


0212-4-7'^ 


128.93 


11.96 


0224+671 


4C +67.05 


132.12 


6.23 


0355+508 


4C +50.11 


150.38 


-1.60 


0415+379 


3C111 


161.68 


-8.82 


0727-115 


PKS 0727-11 


227.77 


3.14 


1730-130 


[HB89] 1730-130 


12.03 


10.81 


2013+370 


[CCS92] 2013+370 


74.87 


1.22 


2037+511 


[HB89] 2037+511 


88.81 


6.04 


2200+420 


[HB89] 2200+420 


92.59 


-10.44 



Note: 'NED Name' gives names under which the sources are 
recorded in the NASA/IPAC Extragalactic Database. 

of AM /M = (-4 ± 6) x 1CT 5 but as we show above this 
comparison actually constrains g P a 2 , not mass. Furthermore 
they used overestimates of both the value and error. They 
used the Wiklind & Combes (1994) measurement which had 
the CO line offset from the HI velocity by —11 kms" 1 ; a cor- 
rected CO measurement (Combes & Wiklind 1996) shows 
there is no measurable offset. Furthermore Varshalovich & 
Potekhin (1996) used the width of the HI line for the mea- 
surement uncertainty. Even allowing for systematic errors 
the true uncertainty is at least a factor of 10 smaller so 
these data in fact establish a limit of order 10 -5 or better. 
This potential for improved limits has prompted the present 
investigation: previous upper limits on change in a are of 
order Aa/a ~ 10 -4 (Cowie & Songalia 1995; Varshalovich, 
Panchuk & Ivanchik 1996). 

We consider three sources for which combined HI 21 cm 
and molecular absorption detections have now been pub- 
lished. Before examining them in detail we look at some ex- 
amples of local Galactic absorption in Section ^ to show that 
the same clouds are detected in both atomic and molecular 
absorption. In Section ^| we summarize the published mea- 
surements of the sources and reanalyze the HI data with 
smaller errors than Varshalovich & Potekhin (1996) to ob- 
tain improved upper limits on the variation of g p o? . The 
interpretations of the results, a comparison with terrestrial 
limits, and future prospects are described in Section H. 



2 SOURCES WITH GALACTIC ABSORPTION 

One possible problem with this approach is that the HI 
and molecular absorptions might result from different clouds 
along the respective lines of sight. We need to show that if 
we take the nearest HI system closest in velocity to a given 
molecular absorption they are associated and not just a ran- 
dom superposition of different absorbing clouds. We inves- 
tigated this by looking at a sample of extragalactic mm- 
wave continuum sources (Liszt & Lucas 1996). Nine of these 
sources have both HCO + absorption measured by Liszt & 
Lucas and HI absorption measured by Dickey et al. (1983). 
We note that these sources are all at low Galactic latitude 
and will therefore experience more absorption than we would 
expect from extragalactic sources intersecting randomly ori- 
ented external galaxies. 

The nine Liszt & Lucas sources with both HI and HCO + 
absorption are listed in Table hi and we plot their HI spectra 
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Table 2. New HI absorption line fits 
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Figure 1. HI and HCO+ spectra of sources with Galactic absorp- 
tion. The HI data (solid lines) are from Dickey et al. (1983) and 
the HCO+ data (dotted lines) are from Liszt &; Lucas (1996). 
The velocity scale is with respect to the rest frequency of the 
respective lines in the local standard of rest frame. 



with the HCO + overlayed in Fig. ^ using data kindly pro- 
vided by the respective authors. They are plotted on com- 
mon velocity scales with respect to the rest frequency f re st 
of the corresponding transition at the local standard of rest 
frame. The velocities v are approximated from the observed 
frequency fobs by v/c = 1 — fobs/ frest where c is the speed 
of light. 

The velocity structure of the HI absorption is complex 
in most of the sources, consistent with the low Galactic lat- 
itudes of these lines of sight. The tendency for most absorp- 
tion to be at negative velocities (as illustrated in Fig. |l|) is a 
consequence of Galactic rotation and the Galactic longitudes 
of the sources. 

Although these lines of sight display many more HI ab- 
sorption systems than for HCO + , our aim is just to test 
the hypothesis that every molecular absorption can be di- 
rectly associated with HI at the closest velocity: to do this 
we fitted velocity components to the HI and molecular data 
independently and the compared the velocity differences sta- 
tistically. 

The procedure adopted for fitting was as follows. As 
described above, we treated the HI and HCO + data inde- 
pendently, so we did not require the same numbers of com- 
ponents in both HI and HCO + systems. We used the VP- 
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z obs 


offset 


1 1 ^ 






( km s ) 




02 18+357 


.684680 


.684495+. 000010 


—32.9+1.8 


1 






.684568+. 000011 


-20.0+1.9 


2 






.684684+. 000006 


0.8+1.0 


3* 


1413+135 


.246710 2 


.246710+. 000004 


0.0+0.9 


1* 






.246770+. 000037 


14.5+8.8 


2 


1504+377 


.671500 3 


.673015+. 000014 


271.6+2.5 


1 






.673236+. 000003 


311.2+0.5 


2* 


U 




.673423+. 000005 


344.7+0.8 


3* 


Li 




.673569+. 000068 


370.8+12 


1 



Notes: z i 3 is the fitted redshift and uncertainty including both 
fitting and velocity scale (0.3 kms -1 ) errors; offset is the velocity 
of the fitted line with respect to the reference redshift z re f taken 
from the references listed in Table ^. ID is the number of the 
component labeled in Fig. ^ with an * to indicate the strongest 
lines used to compare with the molecular data. 

FIT software (Webb 1987) to fit multiple components to 
both the HI and HCO + absorption systems. VPFIT uses an 
unconstrained non-linear least-squares optimization method 
to fit Gaussian/Voigt profiles to the absorption lines in each 
system. We adopted the standard approach of fitting the 
minimum number of components required to obtain a sta- 
tistically acceptable fit, i.e. that the normalized x 2 was close 
to unity. The fitting process was also constrained not to 
introduce any components narrower than the instrumental 
resolution. 

The results of the fitting process are summarized in 
Fig. ^, a histogram of the (HI — HCO + ) velocity differences 
for all systems (solid line) and the closest matching pairs 
(dashed line). The distribution of all pairs shows a peak at 
zero; it is also slewed to negative velocities because HI is 
detected out to larger Galactocentric distances than HCO + . 

To test our hypothesis of a physical association between 
each molecular system and the HI component closest to it in 
velocity, we show in Fig. [] the histogram of closest (nearest- 
neighbour) velocity differences. This distribution is centred 
close to zero (mean Av — 0.4 km s -1 ) and is very narrow 
with a Gaussian dispersion of only 1.2 kms -1 . This is equal 
to the spectral resolution of the HI data (0.6-1.3 kms -1 ) 
(Liszt & Lucas 1996) . To show that the narrow distribution 
of the nearest velocities is significant and does not result 
from choosing nearest pairs from uncorrelated spectra we 
also plot in Fig. |2| the distribution of nearest lines (dotted 
line) where, for each source, we replaced all the HI veloc- 
ities by a random selection of HI lines from all the other 
sources. This shows no peak at zero so we conclude that the 
absorption is from the same clouds to within our velocity 
resolution. 



3 HIGH-REDSHIFT ABSORPTION SYSTEMS 



3.1 Sample of redshifted absorbers 

We found three sources with common extragalactic molec- 
ular and HI absorption detections published: 0218+357, 
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Figure 2. Histogram of the HI — HCO + velocity differences for 
all the sources with Galactic absorption listed above. For each 
source all line pair differences were used to give the solid his- 
togram. The dashed line indicates the nearest HI to each HCO + 
line. The dotted line is the result of a random simulation of the 
nearest pair data to show that the peak for the nearest lines is 
real. 



Figure 3. Spectra of three quasars with high-redshift HI absorp- 
tion. The data are from references 2, 4 and 6 listed in Table pi 
The velocity zero points are taken from the same references. The 
spectra were first normalized to a continuum of unity and then de- 
convolved by fitting multiple absorption components as described 
in the text: the final fit is plotted over the data as a continuous 
smooth line with the centres of each component marked. 



1413+135 and 1504+377. The details of these sources are 
listed below in Table ^[ The sources include 1413+135 which 
was previously discussed by Varshalovich & Potekhin (1996). 
The published HI data are also included in Table [| The typ- 
ical HI error of Az/(l + z) = 1 x 10" 5 (3 kms -1 ) is smaller 
than the value used by Varshalovich & Potekhin (1996), but 
still the dominant term. We show in the next section how we 
improved the HI precision by refitting HI lines to the data 
kindly provided by Dr. C. Carilli. We then compared these 
to the published molecular data to put new limits on the 
variation of g p a 2 . 

3.2 HI absorption systems 

In this section we describe how we obtained more precise 
redshifts for the HI absorption data. A good rule of thumb 
is that the uncertainty in fitting a line position is about one 
tenth of a resolution element. In the case of the HI data for 
1504+377 the resolution was 10 kHz per channel at 848.8 
MHz or about 3 kms -1 , about the error quoted by Carilli 
et al. (1997). Applying the one tenth of a channel criterion 
we might expect to better than this, although we must also 
consider any systematic errors. Checks made by Carilli et 
al. (1997) indicate that the systematic errors in the velocity 
scale are about 0.8 kHz or 0.3 kms" 1 . 

The HI data for three sources are shown plotted in 
Fig. |^. We used to the VPFIT software described above 
to decompose the absorption systems, following the same 
procedure as before, fitting a number of components to each 
system. The resulting fits are plotted over the data in Fig. 



All the components we fitted to each source are listed in Ta- 
ble |^, both in terms of observed redshift and a velocity offset 
with respect to the nominal redshift. The fitting routine gave 
uncertainties less than lkms" 1 for the strong lines; these 
were added in quadrature with the velocity scale uncertainty 
of 0.3 kms" 1 to obtain the errors given in Table |^. 

As with the low redshift galactic absorption we took 
the nearest matches between these HI absorption detections 
and the molecular absorption redshifts listed in Table ^. 

3.3 Molecular absorption systems 

The redshifts of the molecular absorptions are listed in Ta- 
ble ^j. These are — in principle — much more precise than the 
HI data, with velocity uncertainties of 0.1 kms" 1 . In prac- 
tice, the published values were not quoted to this precision 
because they were not aimed at this kind of work. We have 
taken the best estimates of the molecular redshifts and un- 
certainties from the publications as follows. 

In the case of 0218+357 Wiklind & Combes (1995) work 
with respect to the HI redshift of z = 0.68466 (Carilli, Ru- 
pen & Yanny 1993) but it is clear from their Fig. 3 that 
the centre of the molecular absorption (measured from the 
HCO + system) is 4 ± 0.5 kms" 1 higher corresponding to 
z = 0.684680 ± 0.000006, the value we adopt here. 

The velocity of the CO absorption in 1413+135 was 
originally found to be offset from the HI result by Wiklind 
& Combes (1994) but this was subsequently corrected and 
is equal to the HI value of z=0. 24671 (Combes & Wiklind 
1996). We assume the uncertainty in the molecular redshift 
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Table 3. Comparison of Molecular and Atomic Absorption Data 

Name NED Name Molecular z Molecule 



0.684680 ± 0.000006 1 
0.246710 ± 0.000005 3 
0.673350 ± 0.000005 5 
0.673350 ± 0.000005 5 



HCO+(2-l) 
CO(0-1) 
HCO+(2-l) 
HCO+(2-l) 



0218+357 TXS 0218+357 

1413+135 PKS 1413+135 

1504+377 B3 1504+377 
1504+377 

References: 

1: Wiklind & Combes (1995) modified; see text, 
2: Carilli, Rupen & Yanny (1993), 

3: Wiklind & Combes (1994) corrected by Combes & Wiklind (1996), 
4: Carilli, Perlman & Stocke (1992), 
5: Wiklind & Combes (1996), 
6: Carilli et al. (1997). 



Atomic z (HI) Atomic z (HI, new) ^Az/i 



/1 + Z 



0.68466+0. 00004 2 0.684684 ± 0.000006 5 x 10~ 6 

0.24671+0. 00001 4 0.246710 ± 0.000004 5 x 10" 6 

0.67324+0. 00001 6 0.673236 ± 0.000003 

0.67343+0. 00001 6 0.673423 ± 0.000005 



is half the least significant digit, although the potential value 
is much smaller. 

The redshift of the molecular absorption in 1504+377 
quoted by Wiklind & Combes (1996) (z=0.67335) is not con- 
sistent with the HI detections of (Carilli et al. 1997). Carilli 
et al. interpret this difference as either a shift of 15 kms -1 in 
one of the velocity scales or intrinsic differences in the atomic 
and molecular absorption. Their HI redshift has been con- 
firmed by subsequent independent observations (Carilli et al. 
1998, Ap.J., submitted). Our own deconvolution of the HI 
data gives a 4-line system very similar to, but displaced by a 
constant shift from, the HCO + absorption data in Fig. 2(a) 
of Wiklind & Combes (1996). An offset of a single line could 
be interpreted as a case of different clouds causing the ab- 
sorption, but a systematic offset of 4 lines is much more 
likely due to a velocity error in one of the measurements. 
We have not used 1504+377 in our analysis below because 
of this disagreement. 

3.4 Comparison of HI and molecular systems 

Our new more accurate redshift estimates are listed with the 
molecular results in Table ||. In the two sources 0218+357 
and 1413+135 the redshifts agree to within our errors. 
(As discussed above we are not analyzing 1504+377.) We 
can therefore combine the uncertainties in quadrature to 
give 1-sigma upper limits on the redshift differences. These 
are shown in the Table and are |Az/l + z| < 5 x 10~ 6 
(1.5 kms -1 ) for both sources. 

We must still consider the possibility that the molecular 
and atomic absorption arises in different gas clouds along 
the line of sight. This could explain any observed difference. 
However there is no measurable difference between the two 
velocities in our data, so we are probably detecting the same 
gas. The alternative would be that there was a change in the 
frequencies but that in both cases it was exactly balanced by 
the random relative velocity of the two gas clouds observed. 
We consider this very unlikely because of the small 1 km s _1 
dispersion within single clouds shown in Section bl 



4 DISCUSSION 

We can now use the limits to Az/1 + z with the relationship 
given in Section [l] to derive 1-sigma limits on any change in 
y = g p a 2 : \&y/y\ < 5 x 10" 6 at both z = 0.25 and z = 0.68. 



These are significantly lower than the previous best limit of 
1 x 10 -4 by Varshalovich & Potekhin (1996) (it was quoted 
as a limit on nucleon mass, but it actually refers to g p a 2 ). 

As there are no theoretical grounds to expect that the 
changes in g p and a 2 are inversely proportional, we obtain 
independent rate-of-change limits of \g P /g P \ < 2 x 10~ 15 y _1 
and \a/a\ < 1 x 10~ 15 y _1 at z = 0.25 and \g P /g P \ < 
1 x 10 _15 y _1 and \a/a\ < 5 x lO^y -1 at z = 0.68 (for 
Hq = 75kms~ 1 Mpc" 1 as adopted throughout this section). 
These new limits are much lower than the previous 1 sigma 
limit of \a/a\ < 8 x 10~ 15 y" 1 at z ^ 3 (Varshalovich et al. 
1996). 

The most stringent laboratory bound on the time vari- 
ation of a comes from a comparison of hyperfine transi- 
tions in Hydrogen and Mercury atoms (Prestage, Tjoelker 
& Maleki 1995), \a/a\ < 3.7x 10~ 14 y" 1 , and is significantly 
weaker than our astronomical limit. The other strong terres- 
trial limit that we have on time variation in a comes from 
the analysis of the Oklo natural reactor at the present site of 
an open-pit Uranium mine in Gabon, West Africa. A distinc- 
tive thermal neutron capture resonance must have been in 
place 1.8 billion years ago when a combination of fortuitous 
geological conditions enriched the subterranean Uranium- 
235 and water concentrations to levels that enabled spon- 
taneous nuclear chain reactions to occur (Maurette 1972). 
Shlyakhter (1976, 1983) used this evidence to conclude that 
the neutron resonance could not have shifted from its present 
specification by more than 5 x 10" 4 eV over the last 1.8 
billion years and, assuming a simple model for the depen- 
dence of this energy level on coupling constants like a, de- 
rived a limit in the range of \a/a\ < (0.5-1.0) x 10 _17 y _1 . 
The chain of reasoning leading to this very strong bound is 
long, and involves many assumptions about the local con- 
ditions at the time when the natural reactor ran, together 
with modeling of the effects of any variations in electro- 
magnetic, weak, and strong couplings. Recently, Damour & 
Dyson (1997) have provided a detailed reanalysis in order to 
place this limit on a more secure foundation. They weaken 
Shlyakhter's limits slightly but give a 95% confidence limit 
of -6.7 x 10" 17 y _1 < cV/q < 5.0 x 10~ 17 y _1 . However, 
if there exist simultaneous variations in the electron-proton 
mass ratio this limit can be weakened. 

The Oklo results provide stronger limits on the time 
variation of a than the astronomical data; however, the 
astronomical limits have the distinct advantage of resting 
upon a very short chain of theoretical deduction and are 
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more closely linked to repeatable precision measurements of 
a simple environment. The Oklo environment is sufficiently 
complex for significant uncertainties to remain. 

Unlike the Oklo limits, the astronomical limits also al- 
low us to derive upper limits on any spatial variation in a. 
Spatial variation is expected from the theoretical result that 
the values of the constants would depend on local conditions 
and that they would therefore vary in both time and space 
(Damour & Polyakov 1994). The two sources for which we 
derived limits, 0218+357 and 1413+135, are separated by 
131 degrees on the sky, so together with the terrestrial result, 
we find the same values of a to within \Aa/a\ < 3 x 10 -6 in 
three distinct regions of the universe separated by comov- 
ing separations up to 3000 Mpc. Limits on spatial variation 
of g p a m e /m p were previously discussed by by Pagel (1977, 
1983) and Tubbs & Wolfe (1980). We have improved on their 
limits by some 2 orders of magnitude but as our sources are 
at lower redshift, they are not causally disjoint from each 
other. 

With a careful reanalysis of the available data we have 
made a significant improvement on the existing upper lim- 
its on any change of the product of two physical constants 
g p a 2 with cosmic time. The high-redshift measurements are 
now approaching the best terrestrial measurements based on 
the Oklo data. These could be further improved by a fac- 
tor of 2-5 with additional observations that would not be 
difficult to perform such as fitting the atomic and molecu- 
lar data simultaneously, remeasuring the HI absorptions at 
higher spectral resolution, and checking the discrepancy in 
the observations of 1504+377. 
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